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ABSTRACT 

The solitonic sector of the Nambu - Jona-Lasinio model with baryon number 
one is solved in the presence of an external medium. The calculations fully include 
the polarization of both the Dirac sea and the medium as well as the Pauli blocking 
effect. We found that with an increasing density the medium destabilizes the 
soliton. At finite medium density the soliton mass gets reduced whereas the mean 
square baryon radius shows an increase - a swelling of the soliton. At some critical 
density of about two times nuclear matter density there is no localized solution - 
the soliton disappears. 
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1. Introduction 



The chiral symmetry is believed to be a dominant feature of the low-energy 
sector of QCD and in particular, in the structure of the low-lying hadrons. In 
the present absence of a solution to the non-perturbative QCD, different effective 
models have been employed to study the meson and baryon structure. Among 
them the Nambu - Jona-Lasinio model 1 ) seems to play an essential role. From 
one hand the model accounts for the the chiral symmetry breaking - describes a 
structured vacuum and the mesons as elementary excitations of the non-vanishing 
< qq > condensate and from the other - it incorporates the general accepted 
phenomenological picture of the baryon as a bound state of N c valence quarks 
coupled to the quarks of the polarized Dirac sea. Quite encouraging, the model 
shows a significant success in describing the meson sector 2-6 ) as well as the nucleon 
as a baryon number one soliton 7-14 ). 

As it is suggested by the lattice gauge QCD calculations 15 ), at some finite tem- 
perature and/or density one generally expects a restoration of the chiral symmetry 
and a deconfinement, and hence a change of the hadron structure. Following this 
idea to consider the restoration of the chiral symmetry in a medium as a rele- 
vant mechanism for a scale change which modifies the hadron structure, the NJL 
model has been also applied to study medium modifications of the meson proper- 
ties 16-1 ^). Indeed, adding a medium (partially occupied positive continuum) the 
model allows for the (partial) restoration of the chiral symmetry at finite density 
and/or temperature and thus, the model seems to provide an appropriate working 
scheme to study modifications of the meson properties in medium. Therefore, it is 
reasonable to apply the NJL model to study the structure of the more complicated 
objects like baryons in medium since the model allows one to treat exactly the po- 
larizations of both the medium and the Dirac sea as well as the Pauli blocking 
effect. It is the purpose of present work to solve directly the B=l soliton sector of 
the NJL model in a medium of finite density. 
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2. NJL model with a chemical potential 

We start with the lagrangean of the NJL model with scalar and pseudoscalar 
quark-quark couplings 1 ): 

— G — — 
£ = m^dpV - m ** + -[(**) 2 + (*7 5 f*) 2 ]. (2.1) 

Here * describes a quark field of 5 , L r (2)-flavour (u and d quarks), N c = 3 colours 
and average current mass mo = (m u + m^)/2. Introducing auxiliary sigma and 
pion fields by a = — (7**/A 2 and n = — (7*275 r*/ A 2 lagrangean (2.1) can be 
rewritten in the form: 

C = *[z 7 % - g(a + z 75 ^r)]* + —(a 2 + t? 2 ) - (2.2) 

Assuming the boson fields to be classical (zero boson and one fermion loop approx- 
imation) , one can use path integral techniques in euclidean metrics to evaluate the 
corresponding effective action S e j j as a sum of a quark part 

S q eff = Sp\og(J^ + h-^ (2.3) 

and a meson one 

STff = Jd'x^^ + ^-mlUay (2.4) 

Here fi is the chemical potential and h is the single particle Hamiltonian 

a.V 

h = — + Pg (<r + iitf.r ), (2.5) 
z 

5p means the functional and matrix (spin, isospin and colour) trace. In the meson 
part, using the PCAC, the current mass mo is eliminated in favor of the pion mass 
m-jr and the pion decay constant f n . The coupling constant G is related to A 2 by 
G = <7 2 /A 2 where an additional coupling constant g is introduced for convenience. 
Actually, in order to fix our parameter by reproducing the pion decay constant f w , 
g will be fixed later on to the physical coupling constant in order to identify the 
introduced auxiliary fields with the physical pion and sigma mesons. 

Due to the local four-fermion interaction the lagrangean (2.1) is not renor- 
malizable and a regularization procedure with an appropriate cutoff A is needed 
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to make the effective action finite. Actually, only a part of the effective action, 
Sgff([i = 0), coming from the Dirac sea is divergent. Using the proper-time 
regulator 20 ): 



log«^-y A _ 2 ^e-™, (2.6) 

we regularize it. The rest is the contribution coming from the valence and medium 
quarks and it is finite. One can express the energy of the system (see for more 
details refs. 10 ' 18 ) as a sum of a quark part 

E q = Vval e val + Yl e a + TrE%2(^ A ) (2-7) 

M<e a <fi " 

with the proper-time regularization function 

R a (e,A) = — J ^ eX p(-e 2 r), (2-8) 

2 
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and the meson part which up to a trivial factor is given by eq.(2.4). Here M is 
the constituent quark mass, the energies e a are the eigenvalues of the hamiltonian 
(2.5) 

h(f> a = e a <Pa, (2.9) 

and 

r\ va ] = 1 for e va i > 

(2.10) 

Vval = for e vat < 
The baryon density g(r) can be split in valence, medium and sea parts: 

Qval( r ) = Vval<t>lai{*)<t>val{ r ), (2.11a) 

Qmed(r)= E 4(r)0a(r), (2.116) 

M<e a <n 



Qsea(r) = ]- 0U r )Mr)sign(-€ a ), (2.11c) 
z a 

Using the spectral representation of the energies e a the equations of motions 
for the meson fields 

= and 6 -^lL = (2.12) 

00" 07T 
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can be written as 
<J l v < 1 

a 



° = X2 Nc {Yl <f ) a ( PaRi/2( e a, A) - Vval^val^val 

A cm 



Yl Ma} + 

M<e a <fi 



^ = Y2 Nc {Yl ^a^^aRl^a, A) - Vval^val^^val 



_ _ , (2.136) 

M<e a </i 



3. Fixing of model parameters 

Apart from the current mass mo (eliminated using PCAC) the NJL model 
contains also as parameters the coupling constant G (or A 2 in eq.(2.2)) and cutoff A. 
We fix these parameters in the meson (translationally invariant field configuration) 
sector of the vacuum. 

First, in order to express the coupling constant A 2 as a function of the con- 
stituent mass M = g < a > v , where < a > v is the expectation value of the a field 
in vacuum, we use the equation of motion of the sigma field eq.(2.13a) in vacuum 
(the gap equation): 

A 2 = N c g 2 -^R 2 (M, A) + m\ (3.1) 

The equation of motion of the pioneq field eq.(2.136) gives a trivial solution in 
vacuum < jr > v = 0. 

Second, the cut-off A is fixed by reproducing the experimental values of the 
pion decay constant f n and the pion mass in the meson sector of the vacuum. 
This procedure originates from Eguchi 21 ) and is described in detail in ref. 11 ^). This 
fixes the expectation value of the a in vacuum 

< o > v = U (3.2) 

and leads to the following relation between the constituent mass M and the cutoff 
A: 

N c g 2 -\ J ,R 1 (M,A) = 1. (3.3) 
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It is easy to check that this condition insures (in soft pion limit ~ 0) that the 
auxiliary fields a and tx are the physical meson fields. Thus, in this scheme the 
vacuum value of the constituent mass M = gf n is the only free parameter. In 

00 \ 

principle, the empirical values ) of the quark condensate and the quark current 
mass can be used to fix M but actually they still leave a broad range for M. 

At finite chemical potential we use unchanged the values of A 2 and A. We also 
keep g fixed as a particular scale. It should be noted, however, that the solitonic 
solution does not depend on the particular way in which the auxiliary boson fields 
are scaled. 

4. Numerical procedure 

Similar to other works^' 10 ' 12 ' 13 ) we also use the finite quasi-discrete basis and 
numerical method of Ripka and Kahana 23 ) for solving the eigenvalue problem 
(eq.(2.9)) by putting the system in a spherical box of a large radius D. The basis 
is made discrete by imposing a boundary condition at D. It is also made finite by 
restricting the momenta of the basis state to be smaller than a numerical cutoff 
Kmax- Obviously both quantities have no any physical meaning and the results 
should not depend on them. The typical values, which we use, are D ks 20/M and 
K m ax ~ 7M. In addition all checks 10 ' 13 ) concerning the numerical stability of the 
solution with respect to varying the box size and the choice of the numerical cut-off 
have been done such that the actual calculation is completely under control. 

5. Meson sector 

As a first step we look for a translationally invariant solution (meson sector) 
d^a = d^n = at vanishing as well as at finite chemical potential using the 
finite quasi-discrete basis 23 ). Actually, for the solitonic solution later we need the 
translationally invariant meson field configuration, which makes the effective action 
S e jf stationary. To this end we solve the equations of motion eqs. (2. 13a, 2.136) 
together with Dirac equation (2.9) in the case of constant fields. For the pion 
field we get the trivial solution < n >^= (expectation value at finite chemical 
potential) whereas the equation of motion of the sigma field gives the medium- 
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modified value of the constituent quark mass-f M 



= 9 < o- > /i : 




Q 



f + 

M*<e Q <M J 



Mm! 



A 2 



(5.1) 



For a particular value of the chemical potential the baryon density of the 
medium is given by eq.(2.116). In this case the contribution of the Dirac sea 
to the baryon density is exactly zero. It should be also noted that due to the 
boundary conditions on the basis functions, the medium density deviates from a 
constant value in the vicinity of the box surface. In order to diminish this finite 
size effects we always pay attention to have a box size large enough with respect to 
the size of the soliton. To be particular, at fixed medium density the results show 
within 5 % no dependence on both the box radius and the corresponding number 
of the quarks in the box. The stability of our translationally invariant solution 
with respect to finite size effects is illustrated in Fig. 1 where the constituent mass 
M*, calculated in the Kahana-Ripka basis as a function of the baryon density, 
is compared with the results 18 ) obtained using a plane wave basis (no finite size 
effects). As can be seen that up to one and a half nuclear matter density (g n m) 
the two solutions almost coincide and start to deviate slightly at higher density 
values. It shows that the finite basis provides a satisfactorily good description of 
the positive energy continuum (medium). Needless to say that the reduction of the 
constituent mass, found at finite density, indicates a trend towards a restoration 
of the chiral symmetry in medium. 



■t-the asterisk means medium-modified value 
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6. B=l solitonic (nucleon) sector 

As a next step we solve the solitonic sector in medium. We look for a localized 
bound solution of three valence quarks interacting with the Fermi and Dirac sea 
quarks all of a constituent mass M*. We use a numerical self-consistent iterative 
procedure 10 ) based on a method proposed by Ripka and Kahana 2 ^). The proce- 
dure consists in solving in an iterative way the Dirac equation (2.5) together with 
the equations of motion of the meson fields - eqs. (2. 13a, 2.136). The meson fields 
are assumed to be in a hedgehog form 

cr(r) = <r(r) and 7?(r) = ?7r(r) (6-1) 

and are restricted on the chiral circle 

a 2 + 7? 2 = M* 2 /g 2 . (6.2) 

The techniques for the numerical procedure are well known 8 ' 10 ' 12 ) for the case 
of vanishing chemical potential where there is no contribution coming from the 
positive continuum. They can be easily adopted to the case of finite chemical 
potential. As can be seen the latter leads to an additional source (Fermi sea 
quarks) for the meson fields in eqs.(2.13a, 2.136). 

Because of the presence of the soliton both the Fermi sea and the Dirac sea 
get polarized. Subtracting the energy of the unperturbed Fermi and Dirac sea 
(translationally invariant solution) the B=l soliton energy (effective soliton mass) 
is given by the sum of the energy of the valence quarks and the contributions due 
to the polarization of both continua: 

E = Vval^val+ Yl ( e a- e a) + -« £ S{-R3/2( e a» A )-- R 3/2( e a» A )} 

M<e a <^ Z a ' (6.3) 

+ / 7r m 7r 2 J d 3 r{M* /g-a(r)}. 

The energies e° correspond to the translationally invariant solution. 

The density dependence of the B=l soliton mass calculated for two different 
values of the constituent mass M = 420 and M — 465 MeV is shown in Fig. 2. The 
two curves almost coincide and both show a strong reduction of the soliton mass 
at finite density. Hence one should conclude that at least for the constituent mass 
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between 400 and 470 MeV this medium effect seems to be practically independent 
on the particular value of M (vacuum value). For comparison we also present 
the mass of the soliton of the medium-modified a-model lagrangean 18 ). As can 
be seen the latter shows a weaker density dependence. For instance, at g nm it 
was around 20 % whereas the present calculations give about 30%. This number 

0/1 \ 

is close to the estimate coming from the a — u model ) whereas the many-body 
nuclear theory 25 ) gives a smaller number similar to results of the cr-model. It 
should be noted, however, that these calculations 18 ) suffer from the fact that the 
Fermi sea is not included and the Pauli blocking is not taken into account. 

At about two times nuclear matter density we found no solitonic solution. 
It means that at some critical density the inclusion of the medium destroys the 
condition for the existence of a localized solution in the model. Similar effect is 
found 18 ' 2 ^) in other effective models like the Skyrme model with vector mesons 
and the Gell-Mann - Levy cr-model. It is rather tempting to interpret this effect 
indication of a derealization of the nucleon in a dense medium but such 
a conclusion may go out of the scope of the NJL model itself because of the lack 
of confinement. 

It is also interesting to look at the separate contributions to the soliton energy 
coming from the valence quarks as well as from the Fermi and Dirac sea (see Fig. 3). 
At constituent mass M = 420 MeV in the vacuum as well as in medium about 
a half of the soliton energy comes from the polarized Fermi and Dirac sea. With 
increasing medium density both the valence and the summed sea (Fermi+Dirac) 
contribution show a similar reduction. It should be noted, however, that the 
contribution of the Fermi sea alone is negative (opposite sign) and shows a different 
behaviour. Apart from the sign the contribution of the polarized Dirac sea shows 
a similar density dependence. It means that with increasing density the polarized 
Fermi sea trends towards a destabilization of the soliton whereas the polarized 
Dirac sea stabilizes it. As can be seen the contribution of the Dirac sea takes 
over and the summed sea contribution, which is a result of a cancellation of two 
large numbers, is positive. Apparently, at least in the present model calculations 
the contribution, coming from the polarized Fermi sea, plays an important role 
and can not be neglected. Close to the critical density the contributions of the 
polarized Fermi and Dirac compensate each other and the summed sea contribution 
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tends to vanish whereas the valence energy saturates. Since the localized solution 
exists because of the interplay between the valence and the sea contribution the 
vanishing of the latter leads to a destabilization of the soliton. This picture can 
be understood in the context of the partial restoration of the chiral symmetry, i.e. 
reduction of the chiral order parameter and of the constituent mass. The latter 
means a decrease of the gap of the Dirac spectrum. It leads to a reduction of 
the valence quark energy. We observe also that the Fermi and the Dirac sea are 
getting polarized at the same rate. Thus, the summed sea contribution vanishes 
and the soliton gets destabilized. 

Similar to the energy the baryon density distribution of the soliton g so i(r) can 
be split into valence and sea parts: 

Qsol(r) = Vval4>l a i(r)(l>val(r) + {4>a{r)<l>a{r) -^(r)^(r)} 

M*<e a <» (64) 

+ ^^4>\x{ r )^a{r)sign{-e a ). 
L a 

The wave functions 0^ correspond to the translationally invariant solution. 

It is easy to check thet the baryon number of the solitonic solution is one. 
The effect of the medium on the spatial distribution of the soliton is illustrated in 
Fig. 4 where the baryon density of the soliton is plotted for different values of the 
medium density. It is clear that at finite density the soliton gets more extended 
which supports the idea of swelling of the nucleon in medium. Close to two times 
nuclear matter density the solution shows a clear derealization: the distribution 
is rather close to the uniform one. 

As a measure of the spatial extension of the soliton in medium we also calcu- 
lated the mean square radius of the soliton: 

< r 2 > so i= J d\r 2 Qsol (r). (6.5) 

The results are depicted in Fig. 5. As it is expected from the soliton density 
distributions in medium (shown in Fig. 4) the radius increases with increasing 
medium density. At densities close to the critical one it seems to diverge. Similar 
to the mass the relative change of the radius is larger than the estimates 1 ^) in 
(T-model: at g nm the rms radius is increased by 30% whereas the solitonic sector 
of the a-model gives an enhancement of about 20%. 
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7. Summary 

We solved self-consistently the baryon number one solitonic sector of NJL 
model in order to investigate the structure of the solitonic solution with baryon 
number B = 1 immersed in a quark medium of finite density, medium in a self- 
consistent way. The polarization of both continua, namely the Fermi and the Dirac 
sea, as well as the Pauli blocking effect are treated on the same footing. We found 
that the polarization of the medium plays a role as important as the polarization 
of the Dirac sea and with increasing density the medium polarization destabilizes 
the soliton. The B=l solitonic solution obtained gets more extended in medium 
(swelling). The mass gets reduced whereas the mean square radius increases. 
The relative changes of both quantities are larger than the estimates 1 ^) obtained 
from the soliton sector of the medium-modified Gell-Mann - Levy a-model. At 
increasing medium density the soliton show a clear trend to a derealization and 
at some critical value, less than two times nuclear matter density, it disappears. 

The work has been supported by the Bundesministerium fur Forschung und 
Technologie, Bonn (Int. Biiro), by the KFA Jiilich (COSY - Project) and the 
Deutsche Forschungsgemeinschaft. 
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FIGURE CAPTIONS 



1. Constituent quark mass M* as a function of the medium baryon density. 
At vanishing density it is fixed to M = 420 MeV. The solid line represents 
the present calculations with a finite quasi-discrete basis whereas the 
dashed one - the calculations with plane-wave basis. 

2. Soliton mass as a function of the medium baryon density for two different 
values of the constituent quark mass is compared to the mass of the 
soliton (dash-dotted line) of the medium-modified a-model. 

3. Soliton mass and the contributions coming from the valence quarks and 
from the polarized continua as a function of the medium baryon density 
for the constituent mass M = 420 MeV. The separated contributions 
from the polarization of the Fermi and Dirac sea are shown as well. 

4. Density distribution of the soliton in medium for three different medium 
baryon densities. The constituent mass is M = 420 MeV. 

5. Soliton mean square radius as a function of the medium baryon density 
for two different values of the constituent quark mass. 
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